Thermoelectric Power in the Double Exchange Model 
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Employing the Monte-Carlo method and the exact diagonalization, we have investigated the 
temperature dependence of the thermoelectric power (TEP) for the double exchange model in the 
dilute carrier concentration limit. We have found that the TEP follows the Heikes formula in the 
high temperature regime, whereas, in the intermediate temperature regime, the TEP is suppressed 
by the exchange coupling between itinerant electrons and local spins. In the low temperature regime, 
the TEP exhibits an anomalous peak and dip feature near the magnetic transition temperature Tc 
which can be understood based on the magnetic polaron state. We have also found that the TEP, 
in the presence of the magnetic field, shows the positive magnetothermoelectric power near Tc- 
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The strong coupling between the electrical transport 
and the magnetism appears in the systems which show 
the transition from the paramagnetic (PM) insulator 
to the ferromagnetic (FM) metal. Colossal magne- 
toresistance (CMR) doped manganese oxides ffl, di- 
lute magnetic semiconductors 0, and EuB6 [J, [f| 
are typical examples. Simultaneous FM and metallic 
transitions have been well described by the double ex- 
change model in which the itinerant electrons are cou- 
pled to the local spins through the Hund-type exchange 
interaction. [1, 0, H @ 

In these systems, the thermoelectric power (TEP) as 
well as the electrical conductivity shows an anomalous 
behavior due to a strong correlation between conduc- 
tion electrons and local spins. For example, depending 
on the doping ratio, applied magnetic field, or tempera- 
ture, doped perovskite manganese oxides exhibit intrigu- 
ing TEP behaviors such as the sign change, the appear- 
ance of a large peak in the FM phase, and the magne- 
tothermoelectric effect. [H HO, EEl The TEP would be 
affected by various interactions such as electron-electron, 
electron-phonon, Jahn- Teller interaction, etc. The essen- 
tial feature of the TEP in doped perovskite manganese 
oxides, however, can be the explained by the double ex- 
change mechanism. 

It has been well known that the double exchange model 
in the dilute concentration limit produces the magnetic 
polaron. The magnetic polaron is a composite quasi- 
particle of a charge carrier and the magnetic polariza- 
tion field of local magnetic moments induced by the 
carrier. [3, [ll| The magnetic polaron effect enhances the 
effective mass of the carrier and can localize carriers. 
Thus the magnetic polaron state is featured by FM clus- 
ters embedded in the paramagnetic background and by 
the insulating behavior in the resistivity. [1 €31 ] 

In this study, we have examined the TEP behavior for 



the double exchange model in the dilute concentration 
limit. By using the Kubo's formalism, we have obtained 
the temperature dependent TEP and investigated its be- 
havior with respect to the exchange coupling strength. 
Further, we have demonstrated the held dependence of 
the TEP near the magnetic transition temperature Tc- 

The double exchange model which includes the strong 
exchange interaction between itinerant electrons and lo- 
cal spins has the following Hamiltonian form: 
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where Si and denote the spin operators of the local and 
the itinerant electron at the i-th site, respectively. c ia 
(cj CT ) is the annihilation (creation) operator of the elec- 
tron with the a spin state at the i-th site, t is the hopping 
parameter, and Jh represents the exchange coupling be- 
tween local spins and electrons. For simplicity, we have 
considered, in this study, a two-dimensional {2D) square 
lattice system with Ising-type local spins (S z — ±i). 
Previous studies show that the simulation results of the 
double exchange model do not depend much on the di- 
mensionality and lattice type. [It], E| We have neglected 
the electron-electron interaction. The electron-electron 
interaction can give rise to the important effect on the 
electrical and thermal transport properties. To avoid its 
effect, we have considered a system with very low carrier 
concentration. 

In order to calculate the TEP for the double exchange 
model, we have employed the Kubo's formalism of the 
transport coefficients for a system in the presence of both 
electric field and temperature gradient. [H, [2(| HI| The 
thermoelectric power is given by 
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where e is the electron charge and /x is the chemical po- 
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tential. The transport coefficients M^- are denned by 
Mij = lim \Dij5(ijj) 

+ /^V^(o)>}, 0) 

where N is the size of lattice, D u and D12 are the gen- 
eral stiffnesses of the charge and energy, and J\ and J2 
represent their current operators, respectively. 

Using the Monte-Carlo (MC) method combined with 
the exact diagonalization (ED), we have calculated the 
thermal average of the TEP. For the description of ther- 
mal contribution of local spins, we have used the MC 
method based on the standard Metropolis algorithm [23| 
with the periodic boundary condition. The electronic en- 
ergy and states are obtained through the ED method. [24| 

In the high temperature limit, the TEP for the dou- 
ble exchange model becomes nearly constant and simply 
proportional to the entropy per carrier. It is because the 
first term in Eq. [5] becomes zero at high temperature and 
only the second term contributes to the TEP. Then the 
TEP is expressed by the following Heikes formula: [22| 

S T -> 00) -> — = — = -— In , (4 

\e\T \e\ \e\ n 

where s is the entropy per electron, and n is the elec- 
tron density. This feature is well described in Fig. QJa). 
Regardless of the magnitude of Jh, all the curves of the 
TEP converge to the specific horizontal line. This asymp- 
totic line in Fig. [lja) represents the value of the TEP 
(S = —4.8442 fcs/|e|) for the case of dilute concentra- 
tion n = 0.0156. 

In the intermediate temperature regime (0.5 < T < 
10.0), the TEP exhibits different behaviors depending on 
the strength of Jh- As shown in Fig. [TJa), the value 
of TEP is the largest for Jh = 0.0, and becomes re- 
duced with increasing Jh- It is because the strong ex- 
change coupling between an itinerant electron and local 
spin brings down the spin entropy of conduction electrons 
and inhibits the heat current. 

In the low temperature limit (T < 0.5), the TEP for 
the double exchange model shows interesting behavior for 
finite Jh- The inset in Fig. [TJa) reveals that a peak and 
dip feature occurs in the TEP curve for J H 7^ 0. Note 
that the TEP for Jh = 0.0 converges monotonically to 
zero. Transport properties of a metallic system are 
well described by the semi-classical Boltzmann transport 
theory [26|, which yields the T-linear TEP behavior at 
low temperature. On the other hand, the TEP for an in- 
sulating system is known (2?| to be proportional to 1/T. 
Thus it is inferred from the inset in Fig.QJa) that an insu- 
lating phase appears above the dip position temperature 
Td, but disappears with decreasing temperature below 
T D . 

As mentioned above, the double exchange model in the 
dilute carrier concentration limit is illustrated by the con- 
cept of the magnetic polaron.[16] Below the peak position 
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FIG. 1: (color online), (a) The temperature dependence of 
the thermoelectric power S(T) for different Jh parameters in 
the double exchange model. The thin horizontal solid line 
which crosses S = —4.8442 fcs/|e| represents the value of S 
calculated by the Heikes formula. The inset in (a) presents 
the low temperature behaviors of S(T). (b) The magnetic 
field dependence of S(T) for Jn/t = 10.0. For all cases, n 
(electron density) is fixed by n — 0.0156 (one electron in an 
8x8 lattice). 



temperature Tp, the strong exchange coupling produces 
the composite quasi-particle of an itinerant electron and 
local spins, and accordingly the magnetic polaron of FM 
cluster type is formed. As the clusters are far apart, elec- 
trons are trapped and their mobilities are reduced. As 
the clusters are merged and the FM state begins to ap- 
pear near Tq, the mobilities increase again, [la [H[ That 
is, the insulator to metal transition is caused by the per- 
colation of the magnetic polarons. It is thus evident that 
the peak and dip feature of the TEP in Fig.QJa) is closely 
related to the magnetic polarons which are formed at Tp 
and begin to be percolated at Tn. [29| 

Near Tc, the TEP exhibits the strong magnetic 
field dependence as seen in the magnetic and trans- 
port properties. [30] As show in Fig. [Ub), the magni- 
tude of the TEP diminishes with increasing the mag- 
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netic field and its dip position shifts up in temperature. 
Also notable is that the magnetothermoelectric power, 
AS = S(H) — S(0) is positive in this system. This fea- 
ture is also related to the formation of magnetic polarons. 
Because the thermal fluctuation of the local spins is sup- 
pressed in the presence of the magnetic field, the mag- 
netic polarons are easily formed and merged together. 
The percolation gives rise to the reduction of the electri- 
cal resistivity (the negative magnetoresistance) and the 
decline of the insulating phase. Then the emerging metal- 
lic phase suppresses the peak and dip feature near Tc in 
the TEP curve. Hence the TEP which has a negative 
value moves up with increasing the magnetic field. 

In conclusion, we have investigated the temperature 
dependence of the TEP for the double exchange model 
in the dilute carrier concentration limit. The TEP in the 



high temperature regime follows the Heikes formula in- 
dependent on Jh value. The exchange coupling reduces 
the TEP in the intermediate temperature regime. Near 
Tc, the TEP shows anomalous peak and dip feature and 
the positive magnetothermoelectric power, which are ex- 
plained by the concept of the magnetic polaron. 
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